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Abstract 



The introduction and emergence of agriculture into Eastern North America and Europe proceeded very differently: it 
varied in timing, speed, and mechanism — despite similar woodland environments. To resolve the different subsistence 
paths, I employ the Global Land Use and technology Evolution Simulator, a numerical model for simulating demography, 
innovation, domestication, migration and trade within the geoenvironmental context. 1 demonstrate how Europe receives 
a large package of foreign domesticates and converts rapidly. In contrast, Eastern North American trajectories exhibit 
a gradual transition: hunting-gathering and agropastoralism coexist for a long time, and agriculture is integrated slowly 
into the existing subsistence scheme. 
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In the Levant and along the Chinese Yellow and 
Yangtze rivers, the first green revolution occurred between 
10,000, BC and 8000 BC when the predominant hunting- 
gathering life style was replaced by agriculture and pastor- 
alism as the main subsistence style (Willcox, 2005; Kuijt 
& Goring- Morris, 2002; Londo ct al., 2006). Subsequently, 
this new way of actively procuring food resources took 
hold in almost every society around the globe — at differ- 
ent times, with different speed, and different mechanisms. 
Where agriculture had not arisen independently, or where 
the local domesticates where insufficient to serve as staple 
crops, agriculture was imported from the founder regions 
either by prehistoric trade network or by migrants (e.g., 
Wen et al., 2004; Lemmen et al., 2011). Prime exemplary 
regions for non-independent agriculture are Europe and 
Eastern North America (ENA, Price et al., 2001; Smith, 
1997). 

In most of Europe, the transition to farming occurred 
between 6000 BC and 4000 BC, with a continuous tem- 
poral gradient from the southeast (earliest agriculture in 
Greece and the Balkan, Perles 2001) to the northwest 
(southern Scandinavia, Price 2003, and Britainm Whittle 
& Cummings 2007). The most prominent example for a 
large-scale homogeneous agricultural complex is the Lin- 
earbandkeramik culture of central Europe (Liming, 2005). 
None of the crops used had been domesticated locally 
in Europe, but had been imported from the west Asian 
founder regions; everywhere, the emergence of agriculture 
was associated with the appearance of ceramic artifacts. 
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The Neolithic package which was delivered to Europe is 
comprised of food crops and animals that have their ori- 
gins in the wider region of the Fertile Crescent (Flannery, 
1973), it contains wheat, barley, rye, lentils, peas, cattle, 
sheep, goat, and pigs (Willcox, 2005; Luikart ct al., 2001; 
Edwards et al., 2007; Larson et al., 2007; Zeder, 2008). 

The transition to farming occurred much later in North 
America, between 1000 BC and 1000 AD, a period denoted 
for ENA as the Woodland stage. Local domesticates like 
amaranth and sunflower had been cultivated along a dom- 
inating hunting gathering life style, but gradually corn and 
beans were being introduced from Mexico and became the 
staple crops. The Woodland period is characterized by 
continuous pottery artifacts and the gradual increase of the 
dependence on agriculture (Anderson & Mainfort, 2002). 

The success of growing crops is restricted to suitable 
climatic conditions and sufficient soil and water resources. 
With the exception of the dry plateaus along the Western 
cordillcra, North America and Europe exhibit similar cli- 
matic conditions, vegetation types and topography, both 
with a variety of landscapes and features like large rivers, 
mountain ranges, plains and rolling hills. The natural 
vegetation type on the majority of North American and 
European area is temperate forest (Williams ct al., 2000; 
Thompson & Anderson, 2000; Cheddadi ct al., 1997). 

In this study, I employ the Global Land Use and tech- 
nological Evolution Simulator (GLUES) which resolves 
local innovation, migration and cultural diffusion of traits 
(Wirtz & Lemmen, 2003; Lemmen et al., 2011). With 
GLUES I demonstrate that the model is capable to explain 
the global onset of agriculture for many world regions very 
reasonably; here, I restrict my analysis of the global model 
to North America and Western Eurasia. I elaborate on 
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the differential processes leading to agriculture in Europe 
and the Eastern North American Woodland region, and I 
estimate the pre-Columbian population in the Woodland 
region. 

2. Material and Methods 

The Global Land Use and Technological Evolution 
Simulator (GLUES, Wirtz & Lemmen 2003; Lemmen 
2009; Lemmen & Wirtz 2010; Lemmen et al. 2011) de- 
scribes the evolution of regional sociocultural traits con- 
strained by biogeographical setting. Available resources 
exploited by regional societies are based on regional net 
primary production, derived for the past by applying 
Climber-2-simulated temperature and precipitation anom- 
alies (Claussen et al., 1999) on the HAS A climatological 
data base (Leemans & Cramer, 1991), and a bioclimatic 
limit vegetation model (Licth, 1975). The conceptual 
model is outlined below, for details on the algorithms used 
and the mathematical implementation we refer the reader 
to Wirtz & Lemmen (2003) and Lemmen et al. (2011). 

2.1. Sociocultural model 

The sociocultural realm is described by three charac- 
teristic traits of a population and its demographic density. 
The temporal evolution of each trait follows the gradient of 
increased benefit for success (i.e., growth) of its associated 
population (Dieckmann & Law, 1996; Wirtz & Eckhardt, 
1996). The evolution of the distribution of each trait in 
terms of its moments, not a single realization of the trait 
values is simulated. Below, I give a short description of 
the characteristic traits: 

(1) Technology is a trait which describes the efficiency 
of food procurement and the effectiveness of basic health 
care. In particular, it describes the availability of tools and 
weapons, storage and work organization. Included in this 
trait is writing as a means for administration and cultural 
storage, and transport. (2) The share of farming and herd- 
ing activities (energy or time or manpower) with respect 
to the total food sector. (3) The number of agropastoral 
economies available to a regional population; though I do 
not attribute named plants and animals to each economy, 
as an example, a number of four is obtained when pig 
farming, barley and wheat harvesting, and goat herding 
are present in one population. With this setup I follow 
Shcnnan (2001) by directly coupling population dynamics 
with the process of cultural evolution. 

Model parameters are constrained by requiring that the 
emergence of agropastoral centers in five founding regions 
is successfully simulated with appropriate timing as sug- 
gested by Smith (1997, pl2): these are the Fertile Cres- 
cent (7500 BC), central Mexico (7000 BC), south China 
(8500 BC), north China (7800 BC), south central Andes, 
(7000 BC) and the eastern United States (4500 BC) 1 . For 
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each simulation with random parameter variation, a score 
is calculated considering the spatiotemporal distance to 
these centers. The best scoring parameter set from one 
million simulations was selected for the hindcast study in 
this paper. 

The mathematical model conceptually works as follows: 
Within each biogeographically defined region (average size 
127 • 10 3 km 2 ), a Mesolithic hunter-gatherer population is 
incubated at simulation start. The Mesolithic population 
grows with a speed constrained by sociocultural traits, and 
is limited by natural resources. As technologies become 
more powerful, resources are more efficiently exploited 
and first domesticates are brought under control. Demo- 
graphic density and technologically advanced exploitation 
put pressure on the environment and reduce the usable 
resources. Only when the combination of technologies 
and domesticates outcompetes the hunter-gatherer sub- 
sistence with respect to marginal growth rate benefits, a 
local transition to farming occurs. Typically, the local 
transition is completed in less than 200 years, consistent 
with empirical data (e.g., Zvelebil, 1996). 

Trade and migration are vectors of the spread of farm- 
ing from the founding centers to more distant regions. The 
mechanism of trade and migration relies on the concept of 
influence difference (Renfrew & Level, 1979), which de- 
scribes the tendency and direction of exchange. Migration 
occurs with speed proportional to length of the contact 
zone divided by the distance of two neighbors. Trade ex- 
change is proportional to the difference in the products 
traded, and is assumed to reach one hundred times farther 
than migration. Lemmen et al. (2011) recently investig- 
ated the importance of trade and migration as vectors of 
cultural or demic diffusion followed by local adoption in 
the GLUES model; they found that both processes are 
consistent with the radiocarbon record in Europe. 

For the mathematical model, Lemmen et al. (2011) in- 
troduced the age scale "simulated time BC" (simBC) 
to distinguish between empirically determined age mod- 
els and the model time scale; I adopt this notation 
here. Ideally, simBC should be numerically equal to 
year BC as calendar year. GLUES version 1.1.17 was 
used for this study with the standard 685 world re- 
gion Linearbandkeramik (LBK) setup (Lemmen, 2009) 
and no imposed regional climate fluctuation. GLUES is 
free and open source software and can be obtained from 
http://glues.sourceforge.net/. 

2.2. Biogeograpic and bioclimatic forcing 

The sociotechnological model is embedded in the geo- 
graphic and past vegetation context. In GLUES, vegeta- 
tion is measured by its net primary productivity (NPP); 
net primary productivity is derived from a bioclimatic 
limit model (Licth, 1975) on annual mean temperature 
and annual precipitation. To obtain the vegetation of the 
past, I use results from a transient Holocene Climber-2 
(Claussen et al., 1999) intermediate complexity climate 
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Figure 1: Biogeographic 
context for the sim- 
ulation of sociotech- 
nological evolution 
with GLUES. The 
biogeographic vari- 
ables productivity (top 
panel), agropastoral 
potential (middle panel) 
and number of suitable 
species (bottom panel) 
are derived from recon- 
structed temperature 
and precipitation at 
1000 sim BC 
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Figure 2: Regional timing of the transition to agriculture or pastor- 
alism as the dominant subsistence style, hindcasted with the GLUES 
socio-technological simulator. All colored regions have converted 
to agropastoral subsistence, while white coloring indicates that the 
model does not predict the transition there by 1490 AD. 



model simulation as anomalies on the IIASA climatolo- 
gical data base for mean monthly values of temperature 
and precipitation (Leemans & Cramer, 1991). 

Figure 1 shows available natural resources relevant for 
the simulated transition to agriculture: the highest pro- 
ductivity is found in the southeast of North America, 
intermediate values throughout the eastern and central 
United States and throughout Europe, low values in boreal 
Canada, the Rocky Mountains and the Sahara (top panel 
in Figure 1). The figure shows vegetation productivity at 
1000 sim BC, which is representative for the entire simula- 
tion because changes in the climate data set are so small 
that the simulated productivity does not vary much during 
the Holocene. 

Albeit higher in the southeast of North America, the 
vegetation productivity is in a suitable range (500-1500 
grams per square meter and year) in both the ENA and 
European woodland regions, when it is converted to the re- 
lative agropastoral productivity potential (middle panel in 
Figure 1). With the transfer function defined by Wirtz & 
Lcmmcn (2003, their Figure 1) for the conversion from net 
primary productivity to agropastoral productivity poten- 
tial, this potential is high (above 75%) for all of Western 
Europe and the eastern and central United States. 

The bottom panel in Figure 1 shows the relative num- 
ber of suitable species — which is related to the number of 
different subsistence economies — in Europe and Eastern 
North America. The transfer function from net primary 
productivity to local species diversity has a narrow op- 
timal range of 400 — 600 grams per square meter and year 
(Wirtz & Lemmcn, 2003, their Figure 1), which indic- 
ates the open woodland optimum for agriculturally suit- 
able (e.g., annual grasses, grazing herbivores) animal or 
plant resources. The simulated species diversity is very 
high throughout Europe; it is very high in the north east- 
ern and the central United States, but only moderately 
high in the south eastern United States. To account for 
the expected higher diversity in large areas versus small 
areas, already shown almost two centuries ago by Watson 
(1835), I employ a concave species-area relationship to 



estimate from the local resources the contmentally avail- 
able diversity of suitable species (Connor & McCoy, 2001; 
Wirtz & Lemmen, 2003). On a continental scale, the po- 
tential suitable species diversity in Eurasia is four times 
the potential which is to be expected for North America. 

3. Results 

Farming and pastoralism as major subsistence styles 
emerge in the model simulation almost everywhere on 
the globe by 1490 sim AD. Hunting and gathering re- 
main important in the boreal zone of Europe and North 
America (Figure 2), and on large islands like Iceland 
(not shown) and Ireland. The timing of the transition 
to agropastoralism — detected in the model when more 
than half of the regional population rely on farming or 
herding — is depicted by colors in Figure 2. There is a 
timing offset between Europe and ENA, with most of con- 
tinental Europe converting to agropastoralism during the 
period 7000-4000 sim BC, and most of ENA converting 
between 2500 sim BC and 1490 sim AD. This lag in the sim- 
ulated timing of the onset of agriculture between Europe 
and Eastern North America is a consequence of the dif- 
ferent biogeographic setting of the two subcontinents as 
discussed by Wirtz & Lemmen (2003). 

3.1. European agricultural expansion 

In Europe, the earliest simulated elevated population 
levels appear in Greece, Italy, Anatolia, and the northwest- 
ern corner of the Iberian peninsula around 8000 sim BC 
(Figure 3, time from top left to bottom right). By 
6000 sim BC, this elevated population density of approx- 
imately one person per square kilometer is visible through- 
out Europe. A dense population center with up to three 
people per square kilometer is focussed on northern Greece 
and extends throughout the Balkan and Anatolia. At 
5000 sim BC, population intensifies further in these regions 
up to five people per square kilometer. In Italy and on the 
southern Iberian peninsula population intensifies at this 
time. 

By 4000 sim BC, most of continental Europe exhibits 
population densities of a little above two people per square 
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Figure 3: Population 
density (km -2 ) in 
Europe between 
8000 sim BC and 
lOOOsimBC, simulated 
with GLUES. 
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Figure 5: Trajectories of 
socio-technological de- 
velopment in European 
(left column, panel 
a,c,e) and Eastern 
North American (right 
column, panel b,d,f) 
woodlands simulated 
with GLUES and aver- 
aged over all regions in 
Europe and ENA, re- 
spectively. From top to 
bottom: (a,b) fraction 
of agropastoral subsist- 
ence, (c,d) number of 
diverse economies, and 
(e,f) percentage of used 
economies relative to 
the natural potential. 




Figure 4: Population 
density (km -2 ) in 
North America between 
3000 sim BC (a) and 
1000 sim AD (f), sim- 
ulated with GLUES. 
Note the threefold scale 
difference to European 
density (Figure 3, 
ranging up to 6 km -2 
in Europe and up to 
2 km -2 on the North 
American continent. 



kilometer, with a densely populated region of five to six 
people per square kilometer in Greece, the Balkans, Anato- 
lia and surrounding areas. By 3000 sim BC, Italy and the 
southern Iberian peninsula have experienced population 
boosts up to five people per square kilometer, much of 
central Europe is populated at a little less than three 
people per square kilometer. By 1000 sim BC, all along 
the Mediterranean coast population densities reach five 
to six people per square kilometer, and in most of the 
European woodland areas a population density of almost 
three people per square kilometer is supported. 

3.2. North American population expansion 

In contrast to the development in Europe, the simulated 
emergence of high population density in North America, 
supported in part by agricultural or pastoral practice, oc- 
curs later in time and reaches only half of the maximum 
density attained in central Europe. A first increase of pop- 
ulation from the background level of 0.3 people per square 
kilometer population density is visible at 3000 sim BC in 
the area of coastal Virginia and North Carolina (Figure 
4). The lowest population density is simulated along the 



Rocky Mountains with less then 0.1 people per square kilo- 
meter. By 2000 sim BC, the population density in these 
states reaches 0.8 people per square kilometer, popula- 
tion density is slightly elevated above the background level 
in parts of eastern Texas, Arkansas, Missouri, Tennessee, 
Mississippi, Alabama and Georgia. By 1000 sim BC, the 
population density in this southern Eastern North Amer- 
ican region and in Texas reaches up to one person per 
square kilometer, with the highest density in a region on 
the Mississippi river between Arkansas, Missouri, Ken- 
tucky, and Tennessee. 

From IAD, the local population centers of the south- 
ern ENA join with eastern Texas to form a large region of 
dense population extending from Mexico throughout the 
the southern ENA to the Atlantic coast. Highest popula- 
tion density is simulated in the coastal states of North Car- 
olina and Virginia with 1.1 people per square kilometer. 
By 500 sim AD population in this core area continues to 
increase, while in the coastal states to the south, the mid- 
west states to the west and up the fourtieth parallel north 
of this core area population increases to 0.6 people per 
square kilometer. By 1000 sim A. D population increases 
in the entire ENA region south of the fourtieth parallel up 
to 1.5 people per square kilometer at the Atlantic coast 
and in the center, and up to 0.9 in the remaining part 
of ENA. In the region north of the fourtieth parallel, low 
population densities of up to 0.4 people per square kilo- 
mter south of the Great Lakes (up to 0.35 in Canada) are 
hindcasted by the simulation. 

3.3. Trajectories to woodland farming 

Sociocultural trajectories to woodland farming are con- 
trasted for Europe and the American East in Figure 5. 
This figure shows the development of the two socio-cultural 
variables fraction of agropastoralism (panels a, b) and eco- 
nomic diversity (panels c,d), and the realized fraction of 
the economic potential (panels e,f). For both Europe 
and ENA, a 6500 year interval encompassing the trans- 
ition to agriculture is displayed, 7500-1000 sim BC for 
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Europe, and 5000simBC to 1500simAD for ENA. Most 
of the European transition takes place between 6500 and 
4000simBC; in this 2500 year interval 80% of the trans- 
ition occurs (panel a). In contrast, the American trans- 
ition starts after 4000simBC and is not completed during 
pre-Columbian times, i.e., the duration of the transition is 
more than twice as long as for Europe (panel b). Not only 
at the subcontinental level, but also at the individual sim- 
ulation region level, the transition in the Americas takes 
more then twice as long (around 600 years) as the trans- 
ition in most European regions (around 250 years). 

The middle panels in Figure 5 show the trajectories 
of economic diversity — indicative of the number of do- 
mesticates in use. There is a marked difference between 
Europe (panel c) and Eastern North America (panel d). 
In Europe, the model simulates a gradual average increase 
of this quantity, resembling the increase of the importance 
of agriculture. On average, five diverse economies are in 
use by European farmers from 4000simBC, however, this 
number stays approximately constant until the end of the 
simulation. From the same initial value of one half, the 
development in the Americas is much slower: from about 
4000simBC, the number of used diverse economies gradu- 
ally increases but does not exceed 1.5 during the simula- 
tion. The simulated large increase in agropastoral activity 
in ENA relies on only a few different economies. 

While it is clear from this and earlier studies that the 
different maximum value is a consequence of the different 
biogeographic setting, I find a marked difference between 
the two subcontinents Europe and ENA in the relationship 
between actually used economies and the potentially avail- 
able (from local or imported resources) economies; it is 
shown in panels e and f of Figure 5, respectively. The frac- 
tion of used economies increases in Europe from around 
20%-100% before 4000simBC. The low initial fraction in- 
dicates that the abundant resources are hardly used and 
that the system is people-limited: there are not enough 
farmers to put into practice the economic potential dur- 
ing the Neolithic. Only from the Chalcolithic period, the 
variety of resources is efficiently utilized in Europe. In 
North America, already at the beginning of the transition 
to agriculture around 4000 sim BC the full spectrum of the 
economic potential of the landscape is utilized. Here, the 
system is resource-limited, and the import of new econom- 
ies drives the diversification of agricultural practice. 

3. 4- Woodland population development 

Average European population density (panel a) in- 
creased rapidly between 7000 sim BC and 4000 sim BC 
from 0.1 to 3 people per square kilometer. Between 
4000 sim BC and 1000 sim BC, the development is slower 
with another increase of 1.3 people per square kilometer. 
On a total area of 3.1 million square kilometers, this pop- 
ulation density translates to an effective population size of 
13 million at 1000 sim BC. Individual regions show evid- 
ence of diverging evolutions in different European regions, 
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Figure 6: Trajectories 
of demographic develop- 
ment in European (top, 
panel a) and Eastern 
North American (bot- 
tom, panel b) woodlands 
simulated with GLUES 
and averaged over all 
regions in Europe and 
Eastern North America, 
respectively. Both es- 
timates of population 
density (solid line, left 
axis) and population size 
(broken line, right axis) 
are shown. Note the 
time offset of 2500 a 
between the two regions 
and the differing scales 
of the size estimates. 
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Transition start 
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2.8 


1.0 


versity potential 
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2.5 


0.5 
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Early limitation 


people limited 
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Peak population 


13 million @ 


3.5 million @ 




1000 sim BC 


1490 sim AD 



Table 1: Summary of key characteristics of the simulated transition 
to farming in Europe and Eastern North America. 



with population 'explosions' in some areas of the south- 
east (e.g., Anatolia, Greece, Balkan), moderate popula- 
tion increases in the central areas of Europe and very low 
population density in the northern parts (cmp. Figure 3). 

Pre-Columbian population in the American Woodland 
area is hindcasted by the model at around 3.5 million (Fig- 
ure 6a). Of these, one million existed already in the Ar- 
chaic period; the increase during the Woodland and later 
period is 2.5 million. The transition to agriculture and pas- 
toralism is not completed at the end of the pre-Columbian 
period, when hunting and gathering were relevant on half 
the ENA Woodland area (Figure 6b). The American tra- 
jectory to agriculture is characterized by an accelerating 
increase between 4000 sim BC and 1 sim AD, and a linear 
increase between 1 sim AD and 1490 sim AD. Key features 
of the transition on both subcontinents are summarized in 
Table 1. 

4. Discussion 

4-1- Neolithic package and mixed subsistence 

The instant availability of a whole array of diverse eco- 
nomics in Europe can be attributed to arrival of the large 
Neolithic package, which included, for example, emmer, 
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cinkorn, barley, goats, sheep, and pigs (Willcox, 2005). 
The rapid increase in people in the agropastoral sector is 
an indicator that not only the domesticates have been im- 
ported, but possibly also the people familiar with agropas- 
toral life style. This can be interpreted as an outplacement 
of the local hunting-gathering population and replacement 
by migrating farmers, consistent with the demic diffu- 
sion theories put forward by Ammerman & Cavalli-Sforza 
(1973), Sokal ct al. (1991), or recently Haak et al. (2010). 
Alternatively, the fast conversion to agropastoralism can 
be viewed as fast adoption by resident foragers; this mech- 
anism is supported by model simulations with the GLUES 
model by Lemmen et al. (2011): the authors show that 
the transition to agriculture in Southeastern and Central 
Europe is consistent with both a demic or a cultural diffu- 
sion model, and that in either case local adoption played 
a major role. 

Much differently, domesticates in the Eastern North 
American Woodland were comprised of a a small but di- 
verse array of local and imported crops. Locally available 
were goosefoot, squash and sunflower — all part of the East- 
ern Agricultural Complex (Delcourt et al., 1998). Gradu- 
ally, local horticulture was supplemented with the foreign 
domesticates squash, maize, and beans. These were of- 
ten grown together in an assemblage termed the 'Three 
Sisters', and originated from the Mexican founder centers. 
The slow increase and the long transition duration in all 
socio-technological trajectories in the American Woodland 
point to cultural diffusion as the main mechanism of in- 
terregional exchange, and to slow adoption of new ways of 
subsistence. 

4-2. Global and Pre-Columbian Population 

The primary goal of this paper is not to discuss and 
provide a new estimate of pre-Columbian population in the 
Americas. I find it instructive, however, to use estimates 
of population size published in the literature as corner- 
stones for evaluating the model performance with regards 
to simulating population. The global population estimates 
provided by GLUES are within the — admittedly large — 
range of uncertainty of published studies on global prehis- 
toric population (e.g. Coale, 1974; McEvedy & Jones, 1978; 
Birabcn, 2003). At the subcontinental level, the GLUES 
estimate of ten million Europeans (Figure 5a) aligns per- 
fectly with the regional data compiled by McEvedy & 
Jones (1978) or with the regional compilation by Kaplan 
et al. (2011) for 1000 BC. Thereafter, with iron age techno- 
logy, GLUES may underestimate the technological devel- 
opment and population density increases in Eurasia, be- 
cause the efficiency gains by metal working and increased 
social stratification are not implemented in the cultural 
model yet. For the Americas, Africa, and Australia, the 
simulations are valid until the beginning of colonialism, 
which is not implemented in the model yet. 

For pre-Columbian America, the uncertainty within and 
between population studies is very large. Published pop- 
ulation estimates show more than an order of magnitude 



range between 10 and 150 million on the American contin- 
ent; a value of around 60 million was assumed by Nevle & 
Bird (2008) 's study on the impact of prehistoric popula- 
tion. Of these, probably only 10% lived north of Mexico, 
while most of the population was concentrated in the Aztec 
and Incan empires, and in the Amazon basin (Denevan, 
1992). For the American Woodland area, a rough estim- 
ate of four million individuals seems appropriate based on 
Nevle & Bird (2008) and Denevan (1992). 

Only recently, Peros et al. (2010) calculated postglacial 
and precontact population growth rates for North Amer- 
ica based on the Canadian Archaeological Radiocarbon 
Database (CARD, Morlan 2005). Their calculation is a 
maximum pre-Columbian population consisting of 1.8-4.7 
million individuals, with a best estimate of 2.5 million. 
The GLUES model hindcasts 3.5 million individuals at 
1490simAD for Eastern North America (Figure 5d), for 
the whole of North America, the simulation yields 6.5 mil- 
lion people. Given the uncertainties in obtaining the pop- 
ulation estimates, GLUES provides a realistic population 
estimate for North America. 

4-3. Model critique 

The good representation of pre-Columbian Eastern 
American population seems surprising, given the simplistic 
approach and parameter uncertainties associated with the 
socio-cultural model GLUES. What the model does is to 
deterministically translate the geographical setting and ve- 
getation resources into a potential for cultural evolution, 
with universal parameter settings for all world regions; 
then the cultural developments within each region and the 
interaction by trade and migration between regions pro- 
duce estimates of demography, subsistence and technolo- 
gical advance. I interpret the success of this kind of simu- 
lation in the way that the large scale patterns of cultural 
evolution worldwide are based on very similar functional 
dependencies on biogeographic context and cultural inter- 
actions, despite the wide range of observed differences in 
cultural traits. 

The model could be improved even further if good 
chronologies were available for a multitude of archaeolo- 
gical sites and with concurrent clear interpretation. For 
Europe, Lemmen et al. (2011) used the data compilation 
of 765 radiocarbon dates with Neolithic attribution from 
Pinhasi et al. (2005); they found a model-data mismatch 
of ±500 years on average. For the American woodland, 
I here use a selection of radiocarbon dates from CARD 
(Morlan, 2005): 3705 dates carry the designation 'Wood- 
land'; the central 95 percentile of these dates is confined 
to the period 1000 BC to 1500 AD, with maximum fre- 
quency in the 11th century AD. When I apply this stat- 
istic to the GLUES model regions located within the area 
where radiocarbon dates with Woodland signature have 
been identified, the 23 model regions' timing of agropas- 
toral onset is represented by a broad distribution between 
3500simBC and 1490simAD, with maximum occurrence 
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between 2000simBC and lOOOsimAD. For ENA, there 
seems to be a model bias towards earlier dates. 

Compilations of radiocarbon dates, however, are prob- 
lematic on their own, due to the differences in the dated 
material, the differences in calibration procedure, different 
labs employing diverse quality checks and analysis meth- 
odology, and subjective choice of sites to include in the 
data base; not the least, the calibration itself introduces a 
dating uncertainty (e.g. Stuiver et al., 1998). The greater 
problem to overcome in a model-data comparison as done 
in this study is, however, the mismatch in the definition of 
the transition to agriculture. In the simulation, a quantit- 
ative measure for the timing of the transition to agropas- 
toralism is available, defined as the time (simBC) when 
more than half of the food production is performed by 
agropastoralists. This measure is not available in the ar- 
chaeological record, where — amongst others — remains of 
crop or domestic animal taxa, agricultural tools, picto- 
grams of agricultural technology give a hint of the exist- 
ence, not the quantity, of agricultural activity. When this 
is associated with a ceramic or stone tool characteristic, 
agricultural activity may be inferred for sites without dir- 
ect evidence of agriculture. There is no measure of farming 
activity from the archeological record which can be directly 
compared to the simulation results. 

Using the regional patterns of the appearance of agro- 
pastoral evidence in the radiocarbon record and the model, 
and using indirect measures like population density or size, 
we can, however, make quantitative comparisons between 
data and models. In doing this, the most interesting find- 
ing should be the failure of a model to represent the obser- 
vations: where the model hindcast deviates from the his- 
toric pathway, that is an indication of the need for explan- 
ation in the cultural realm. One such example could be the 
unusually (with respect to the natural capacity of the en- 
vironment) low population density estimated for hunting- 
gathering societies in North America when compared to 
Europe. This may have led the model-which includes a 
higher population according to the natural capacity — to 
simulate the introduction of agriculture in the ENA Wood- 
land economy much earlier than has been observed (Fig- 
ure 4). Another example of a model shortcoming is the 
incomplete transition to agriculture in North England, Ire- 
land or Scandinavia, which is observed in the radiocarbon 
record (e.g., Price, 2003), but not simulated by GLUES 
(Figure 3). 

The model presented here operates at the large level 
of a society. Necessarily, it does not consider nor ad- 
dress individual cultural choice and agency, i.e. of how 
people actively create and change the society they live in 
(Dornan, 2002). Should the model be required to include 
agency at the collective level of a society, as suggested by 
Shanks & Tilley (1987)? I believe not, because by expli- 
citly excluding agency and simulating only the non-agency 
part of socio-technological evolution, the comparison of 
the model prediction and the historic record reveals where 
agency was important. In other words, the model pre- 



dictions should be interpreted as providing 'a historical 
null hypothesis. Its predictions can be taken as requiring 
no special explanation, and its failures can be taken as 
evidence of rare events that had significant and long-lived 
consequences' (Ackland et al., 2007). 

The current model version lacks transport via sea routes 
and it has no representation of rivers. The planned future 
implementation of sea connections would increase the real- 
ism of simulations on the Mediterranean islands, in Great 
Britain (which is connected to France by an artificial land 
bridge in the model) and Ireland. While we could success- 
fully simulate the Indus Valley transition to agriculture 
(Lemmen & Khan, 2011), the implementation of rivers 
would make the model applicable to a river-based soci- 
ety like Egypt, where the moisture source is far from the 
agricultural area. 

5. Conclusion 

I presented a global numerical model of regional trans- 
itions to agriculture. This model is capable to realistically 
simulate a potential history which can then be compared 
to archaeological evidence. The validity of the simulation 
is exemplified by the model prediction of realistic values for 
population density: ten million Europeans by 1000 sim BC, 
and 6.5 million people on North America (of these 3.5 mil- 
lion in the Woodland area) by 1490 sim A. D. 

Despite the similar environmental setting of the Amer- 
ican and European woodlands, and despite the common 
allochtonous origin of the major staple crops, the trans- 
ition to agriculture and pastoralism proceeded very differ- 
ently in these two regions (Table 1). In Europe, an early 
and rapid adoption of the new lifestyle is simulated with 
the socio-technological simulator GLUES, consistent with 
the archaeological record. A radial expanse from the East- 
ern Mediterranean with a speed comparable to the wave- 
of-advance model is realistically simulated, almost all of 
Europe has adopted agropastoralism by 4000 sim BC. The 
time span of the transition is short and rather people- 
limited than resource limited. 

If the same global parameter set is used to simulate the 
onset of agriculture in the American woodland, the model 
predicts a qualitatively different transition. It is based 
on long-term coexistence of both hunting-gathering and 
agropastoral subsistence, on a mixture of local and foreign 
domesticates, transmitted by cultural diffusion rather than 
demic diffusion. The transition can be characterized as 
resource-limited; it occurs late in time — starting around 
2000 sim BC — and has not completed by 1490. 

Acknowledgments.. This study was partly funded by the 
German National Science Foundation (DFG priority pro- 
ject 1266 Interdynamik) and the PACES programme of 
the Helmholtz Association of German Research Centers. 
Neither of the sponsors had any influence on the content 
or the submission decision of this manuscript. I thank 
K. W. Wirtz for his contributions to and innovations in the 



7 



GLUES model and fruitful comments on prior versions this 
manuscript. I am grateful to D. G. Anderson for his invit- 
ation to present a prior version of this paper at the 75 th 
annual meeting of the Society for American Archaeology 
(St. Louis, April 14-18, 2010) in the session "Trajectories 
to Complexity in Woodland Environments: Eastern North 
America and Temperate Europe compared" . 

References 

Ackland, G. J., Signitzer, M., Stratford, K., & Cohen, M. H. (2007). 
Cultural hitchhiking on the wave of advance of beneficial tech- 
nologies. Proceedings of the National Academy of Sciences, 104, 
8714. 

Ammerman, A. J., & Cavalli-Sforza, L. L. (1973). A population 
model for the diffusion of early farming in Europe. In C. Renfrew 
(Ed.), The explanation of cultural change (pp. 343-57). London: 
Duckworth. 

Anderson, D., & Mainfort, R. (Eds.) (2002). The Woodland South- 
east. University of Alabama Press. 

Biraben, J.-N. (2003). The rising numbers of humankind. Population 
et societes, 394, 1—4. 

Chcddadi, R., Yu, C, Guiot, J., Harrison, S. P., & Prentice, I. C. 
(1997). The climate of Europe 6000 years ago. Climate Dynamics, 
13, 1-9. 

Claussen, M., Kubatzki, C, Brovkin, V., Ganopolski, A., 
Hoelzmann, P., & Pachur, H.-J. (1999). Simulation of an abrupt 
change in Saharan vegetation in the mid-Holoccnc. Geophysical 
Research Letters, 26, 2037-2040. 

Coale, A. J. (1974). The history of the human population. Scientific 
American, 231, 40-51. 

Connor, E. F., & McCoy, E. D. (2001). Species-area relationships. 
Encyclopedia of biodiversity , 5, 397-411. 

Delcourt, P. A., Delcourt, H. R., Ison, C. R., Sharp, W. E., & Gre- 
million, K. J. (1998). Prehistoric human use of fire, the eastern 
agricultural complex, and Appalachian oak-chestnut forests: pa- 
leoecology of Cliff Palace Pond, Kentucky. American Antiquity, 
63, 263-278. 

Denevan, W. M. (1992). The native population of the Americas in 

1492. University of Wisconsin Press. 
Dieckmann, U., & Law, R. (1996). The dynamical theory of coevolu- 

tion: a derivation from stochastic ecological processes. Journal of 

Mathematical Biology, 34 , 579-612. Adap. 
Dornan, J. (2002). Agency and archaeology: Past, present, and 

future directions. Journal of Archaeological Method and Theory, 

9, 303-329. 

Edwards, C. J., Bollongino, R., Schcu, A., Chamberlain, A., Tresset, 
A., Vigne, J.-D., Baird, J. F., Larson, G., Ho, S. Y. W., Heupink, 
T. H., Shapiro, B., Freeann, A. R., Thomas, M. G., Arbogast, 
R., Arndt, B., Bartosiewicz, L., Benecke, N., Budja, M., Chaix, 
L., Choyke, A. M., Coqueugniot, E., Dohle, H.-J., Goldner, H., 
Hartz, S., Helmer, D., Herzig, B., Hongo, H., Mashkour, M., 
Ozdogan, M., Pucher, E., Roth, G., Schade-Lindig, S., Schmolckc, 
U., Schulting, R. J., Stephan, E., Uerpmann, H.-P., Voros, I., 
Voytek, B., Bradley, D. G., & Burger, J. (2007). Mitochondrial 
DNA analysis shows a Near Eastern Neolithic origin for domestic 
cattle and no indication of domestication of European aurochs. 
Proceeding of the Royal Society B, 274, 1377-1385. 

Flannery, K. V. (1973). The origins of agriculture. Annual Review 
of Anthropology, 2, 271-310. 

Haak, W., Balanovsky, O., Sanchez, J. J., Koshel, S., Za- 
porozhchenko, V., Adlcr, C. J., Der Sarkissian, C. S. I., Brandt, 
G., Schwarz, C, Nicklisch, N., Dresely, V., Fritsch, B., Balan- 
ovska, E., Villems, R., Meller, H., Alt, K. W., Cooper, A., & the 
Genographic Consortium (2010). Ancient DNA from european 
early Neolithic farmers reveals their Near Eastern affinities. Pub- 
lic Library of Science Biology, 8, el000536. 

Kaplan, J. O., Krumhardt, K. M., Ellis, E. C, Ruddiman, W. F., 
Lemmen, C, & Klein Goldcwijk, K. (2011). Holocene carbon 



emissions as a result of anthropogenic land cover change. The 
Holocene, 21, 775-791. 

Kuijt, I., & Goring-Morris, N. (2002). Foraging, farming, and social 
complexity in the pre-pottery Neolithic of the southern Levant: A 
review and synthesis. Journal of World Prehistory, 16, 361-440. 

Larson, G., Albarella, U., Dobney, K., Rowley-Conwy, P., Schibler, 
J., Tresset, A., Vigne, J.-D., Edwards, C. J., Schlumbaum, A., 
Dinu, A., Balacsescu, A., Dolman, G., Tagliacozzo, A., Mana- 
seyran, N., Miracle, P., Van Wijngaarden-Bakker, L., Masseti, M., 
Bradley, D. G., & Cooper, A. (2007). Ancient DNA, pig domest- 
ication, and the spread of the neolithic into Europe. Proceedings 
of the National Academy of Sciences, 104, 15276. 

Leemans, R., & Cramer, W. P. (1991). The IIASA database for mean 
monthly values of temperature, precipitation and cloudiness of a 
global terrestrial grid.. Research report International Institute of 
Applied Systems Analyses, Laxenburg. 

Lemmen, C. (2009). World distribution of land cover changes during 
pre- and protohistoric times and estimation of induced carbon 
releases. Geomorphologie: relief, processus, environnement, 4, 
303-312. 

Lemmen, C, Gronenborn, D., & Wirtz, K. W. (2011). A simula- 
tion of the Neolithic transition in Western Eurasia. Journal of 
Archaeological Science, 38, 3459-3470. 

Lemmen, C, & Khan, A. (2011). A simulation of the neolithic trans- 
ition in the indus valley. In L. Giosan (Ed.), Past climates, land- 
scapes and people. American Geophysical Union volume submit- 
ted. 

Lemmen, C, & Wirtz, K. W. (2010). Socio-technological revolu- 
tions and migration waves re-examining early world history with 
a mathematical model. In D. Gronenborn, & J. Petrasch (Eds.), 
The Spread of the Neolithic to Central Europe (pp. 25-38). Mainz, 
Germany: Romisch-Germanisches Zcntralmuscum volume 4 of 
RGZM Tagungen. 

Licth, H. (1975). Modeling the primary productivity of the world. 
Primary Productivity of the Biosphere, 14, 237-263. 

Londo, J. P., Chiang, Y.-C, Hung, K.-H., Chiang, T.-Y., & Schaal, 
B. A. (2006). Phylogcography of Asian wild rice, Oryza rufipogon, 
reveals multiple independent domestications of cultivated rice. 
Oryza sativa Proceedings of the National Academie of Science 
USA, 103, 9578-9583. 

Luikart, G., Gielly, L., Excoffier, L., Vigne, J.-D., Bouvet, J., & 
Tabcrlet, P. (2001). Multiple maternal origins and weak phylogeo- 
graphic structure in domestic goats. Proceedings of the National 
Academy of Sciences, 98, 5927-5932. 

Liining, J. (2005). Bandkeramische Hofplatzc und absolute Chro- 
nologie der Bandkcramik. In J. Liining, C. Friedrich, & A. Zi- 
mmcrmann (Eds.), Die Bandkeramik im 21. Jahrhundert (pp. 
49—74). Rahden (Westf): Symposium in der Abtei Brauweiler bei 
Koln vom 16.9.-19.9. 2002, Internationale Archaologic: Arbcits- 
gemeinschaft, Symposium, Kongrcss volume 7. 

McEvedy, C, & Jones, R. (1978). Atlas 0} world population history. 
Harmondsworth, Middlesex, England: Penguin Books Ltd. 

Morlan, R. E. (2005). Canadian Archaeological Radiocarbon Data- 
base (CARD). Database Canadian Museum of Civilization and 
Canadian Archaolcogical Association. Accessed 2010-04-01. 

Nevle, R. J., & Bird, D. K. (2008). Effects of syn-pandemic fire re- 
duction and reforestation in the tropical Americas on atmospheric 
C02 during European conquest. Paleogeography Palaeoclimatology 
Palaeoecology, 264 , 25-38. 

Perles, C. (2001). The Early Neolithic in Greece: the first farming 
communities in Europe. Cambridge University Press. 

Peros, M. C, Munoz, S. E., Gajewski, K., & Viau, A. E. (2010). Pre- 
historic demography of North America inferred from radiocarbon 
data. Journal of Archaeological Science, 37, 656-664. 

Pinhasi, R., Fort, J., & Ammerman, A. J. (2005). Tracing the origin 
and spread of agriculture in Europe. PloS Biology, 3, 410. 

Price, T. D. (2003). The arrival of agriculture in Europe as seen from 
the north. In The widening harvest: the Neolithic transition in 
Europe: looking back, looking forward (pp. 273-294). Boston, MA, 
USA: Archaeological Institute of America volume 6 of Colloquia 
and conference papers. 



8 



Price, T. D., Bentley, R. A., Liming, J., Groncnborn, D., & Wahl, 
J. (2001). Prehistoric human migration in the Lincarbandkcramik 
of central Europe. Antiquity, 75, 593-603. 

Renfrew, C, & Level, E. (1979). Exploring dominance: predict- 
ing polities from centers. In Transformations: Mathematical Ap- 
proaches to Culture Change (pp. 145-166). New York: Academic 
Press. 

Shanks, M., & Tilley, C. Y. (1987). Re- constructing archaeology: 

Theory and practice. New Studies in Archaeology. Cambridge: 

Cambridge University Press. 
Shcnnan, S. (2001). Demography and cultural innovation: a model 

and its implications for the emergence of modern human culture. 

Cambridge Archaeological Journal, 11, 5-16. 
Smith, B. D. (1997). The initial domestication of Cucurbita pepo in 

the Americas 10,000 years ago. Science, 276, 932-934. 
Sokal, R. R., Oden, N. L., & Wilson, C. (1991). Genetic evidence for 

the spread of agriculture in Europe by demic diffusion. Nature, 

351, 143-145. 

Stuiver, M., Reimer, P. J., Bard, E., Beck, J. W., Burr, G. S., 
Hughcn, K. A., Kromer, B., McCormac, G., van der Plicht, J., 
& Spurk, M. (1998). INTCAL98 radiocarbon age calibration, 
24,000-0 cal BP. Radiocarbon, 40, 1041-1083. 

Thompson, R. S., & Anderson, K. H. (2000). Biomes of western 
North America at 18,000, 6000 and 14 c yr bp reconstructed 
from pollen and packrat midden data. Journal of Biogeography , 
27, 555-584. 

Watson, H. C. (1835). Remarks on the geographical distribution of 
British plants: chiefly in connection with latitude, elevation, and 
climate. London: Longman, Rees, Orme, Brown, Green, and 
Longman. 

Wen, B., Li, H., Lu, D., Song, X., Zhang, F., He, Y., Li, F., Gao, Y., 
Mao, X., Zhang, L., Quian, J., Tan, J., Jin, J., Huang, W., Deka, 
R., Su, B., Chakraborty, R., & Jin, L. (2004). Genetic evidence 
supports demic diffusion of han culture. Nature, 431, 302-305. 

Whittle, A., & Cummings, V. (2007). Going Over: the Mesolithic- 
Neolithic transition in North- West Europe volume 144 of Proceed- 
ings of the British Adademy. Oxford University Press. 

Willcox, G. (2005). The distribution, natural habitats and availab- 
ility of wild cereals in relation to their domestication in the near 
east: multiple events, multiple centres. Vegetation History and 
Archaeobotany , 14, 534—541. 

Williams, J. W., Webb III, T., Richard, P. H., & Newby, P. (2000). 
Late Quaternary biomes of Canada and the eastern United States. 
Journal of Biogeography, 27, 585—607. 

Wirtz, K. W., & Eckhardt, B. (1996). Effective variables in eco- 
system models with an application to phytoplankton succession. 
Ecological Modelling, 92, 33-53. 

Wirtz, K. W., & Lemmcn, C. (2003). A global dynamic model for 
the Neolithic transition. Climatic Change, 59, 333-367. 

Zeder, M. A. (2008). Domestication and early agriculture in the 
Mediterranean basin: Origins, diffusion, and impact. Proceedings 
of the National Academy of Sciences, 105, 11597-11604. 

Zvclcbil, M. (1996). The agricultural frontier and the transition to 
farming in the circum-Baltic region. In D. R. Harris (Ed.), The 
Origins and Spread of Agriculture and Pastoralism in Eurasia 
(pp. 323-345). London: UCL Press. 



9 



